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Abstract 
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^^ 

A study is presented of the process jp — > XY, where there is a large rapidity 
k>( I gap between the systems X and Y. Measurements are made of the differential cross 

^ ' section as a function of the invariant mass M^ of the system produced at the photon 

vertex. Results are presented at centre of mass energies of (W) = 187 GeV and 
(W) = 231 GeV, both where the proton dominantly remains intact and, for the first 
time, where it dissociates. Both the centre of mass energy and the M^ dependence 
of HERA data and those from a fixed target experiment may simultaneously be 
described in a triple-Regge model. The low mass photon dissociation process is 
found to be dominated by diffraction, though a sizable subleading contribution is 
present at larger masses. The pomeron intercept is extracted and found to be 
0^(0) = 1.068 ± 0.016 (stat.) ± 0.022 (syst.) ± 0.041 (model), in good agreement 
with values obtained from total and elastic hadronic and photoproduction cross 
sections. The diffractive contribution to the process jp -^ Xp with M'^/W'^ < 0.05 
is measured to be 22.2 it 0.6 (stat.) ± 2.6 (syst.) ± 1.7 (model) % of the total jp 
cross section at (W) = 187 GeV. 
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1 Introduction 

At the HERA electron-proton collider the bulk of the cross section corresponds to pho- 
toproduction, in which a beam electron is scattered through a very small angle and a 
quasi-real photon interacts with the proton. For such small photon virtualities the domi- 
nant interaction mechanism takes place via the fluctuation of the photon into a hadronic 
state [jl| which interacts with the proton via the strong force. High energy photoproduc- 
tion therefore exhibits broadly similar characteristics to hadron-hadron collisions, but with 
cross sections typically reduced by factors of order the fine structure constant. Whilst the 



regime of asymptotic freedom in strong interactions can be described using perturbative 
QCD, the overwhelming majority of the cross section in which no hard scale is present 
remains far from understood. The study of soft interactions between photons and protons 
provides an excellent environment in which to study both the hadronic manifestation of 
the photon and the dynamics of the long range component of the strong force. 

The dependence on centre of mass energy of elastic and, via the optical theorem, total 
hadron-hadron cross sections has been remarkably well described in a large kinematic 
domain by Regge phenomenology 0. In this framework, interactions take place via the 
t-channel exchange of reggeons related to mesons and of the leading vacuum singu- 
larity, the pomeron (P) [Q. The pomeron is the mediator of diffractiveQ processes. At 
asymptotically large energies, pomeron exchange dominates the elastic channel, such that 
both elastic and total cross sections display a slow increase with centre of mass energy. 
Interactions in which one or both of the hadrons dissociates to higher mass states also 
occur P|. Such processes are characterised by the presence of large regions of rapidity 
space in which no hadrons are produced and are dominated by diffractive exchange at 
large centre of mass energy y/s and small dissociation mass M. The inclusive dissociation 
mass distribution may be treated via Mueller's generalisation of the optical theorem , 
such that an appropriate Regge description involves diagrams that contain three-reggeon 
couplings. 

Experimental results on dissociative processes and their theoretical description are exten- 
sively covered in a number of review articles [3-[TT[| . At sufficiently high energy, differen- 
tial dissociation cross sections da/dM^ are approximately independent of s and fall like 
1/M^ [p!2| , p]3[| . The dependence on the 4-momentum transfer squared t is approximately 
exponential at small t. The highest energy pp experiments |T^,[T^ have confirmed this 
behaviour, but have reported weaker dependences of integrated diffractive cross sections 
on s than is predicted by Regge models based on fits to total and elastic cross sections. 



This is often taken to be the first hint of the infiuence of unitarity bounds |]T6| on diffrac- 
tive cross sections and a number of mechanisms that unitarise the cross section have been 
proposed ^7\, ^ . 



As in the case of hadron-hadron interactions, distinct elastic and dissociative photopro- 
duction subprocesses are usually distinguished. Each of these may be discussed in terms of 
the diagram for the process 7p -^ XY shown in figure |l[ The elastic (EL) case, '-/p -^ Vp, 
corresponds to the situation where the system X consists of a vector meson state, V = p°, 
u, (f), . . ., and y is a proton. For single photon dissociation (GD), 7p — > Xp, the proton 
remains intact and the photon dissociates to a continuum of states with invariant mass 
M^. Single proton dissociation (PD), jp —>■ VY, describes the case where a vector meson 
is produced at the photon vertex, and the proton dissociates to a system of mass M^. 
Finally, in the case of double dissociation (DD), 7p —>■ XY, both interacting particles 
dissociate. 



Results from HERA on the total 'jp cross section [19,20| and exclusive light vector meson 



^In this paper, the word 'diffractive' is used synonymously with i-channel pomeron exchange; s-channel 
approaches to diffraction are discussed in 



photoproduction |pT|-|2^ are well described in the Regge picture in terms of a single 
pomeron trajectory, consistent with that obtained by Donnachie and Landshoff in fits to 
pp and pp total cross sections |2^. Differential cross sections for the GD process have been 



measured in a fixed target photoproduction experiment |^ . In [|T^ overall cross sections 
for the EL and the sum of the GD, PD and DD processes were measured at an average 
centre of mass energy (W) = 180 GeV and in pO], the EL, GD and PD contributions 



were presented as a function of the unmeasured DD cross section at (W) = 200 GeV. 
No measurements have yet been made at HERA that fully unfold the PD, GD and DD 
components or that do not rely on model dependent extrapolations of dissociation mass 
spectra. 

In this paper, the mass M^ is reconstructed directly in the central components of the 
HI detector. In addition to the central detectors, components of HI that are sensitive 
to hadronic activity very close to the outgoing proton direction are used to constrain the 
mass My . This leads to measurements of da/dM^ , both for events in which the proton 
dominantly remains intact (closely corresponding to the EL and GD processes) and where 
it dissociates to low mass states (corresponding to the PD and DD processes). In sec- 
tion]^, the processes measured are defined at the level of final state hadrons and kinematic 
variables are introduced. Section ^ briefly describes the HI detector and section |^ out- 
lines the Monte Carlo simulations used in the measurement. Event selection, kinematic 
reconstruction methods and other details of the experimental procedure are covered in 
section |^. In section ^ the measured cross sections are presented, the GD process is de- 
composed according to a triple-Regge prescription and the dissociation mass and centre 
of mass energy dependences of the diffractive contribution are discussed. 



2 Cross section definitions and kinematics 

Dissociation mass spectra have traditionally been measured by tagging leading final state 
protons and inferring the mass of the dissociation products from the 4- vectors of the beam 
and tagged protons. In this analysis a complementary approach is taken, making use of 
the fact that, where a colourless exchange takes place, there is generally an associated 
region in rapidity space that contains no final state particles (a 'rapidity gap'). The size 
of the gap is related to the masses of the final state systems produced at each vertex. 

Each event is considered in terms of the generic quasi two-body process, jp -^ XY, 
illustrated in figure |T]. The two systems X and Y are separated by the largest rapidity 
gap, with rapidity calculated in the photon-proton centre of mass frame. Y is the system 
closest to the proton beam direction. This scheme of event decomposition provides a 
means of defining hadron level cross sections without assumptions as to the interaction 
mechanism. However, where a colourless exchange does not take place, the location of 
the largest gap is determined by random fluctuations, the average gap size is small and 
at least one of the systems X or Y generally has a large mass. 

The kinematics of the interaction are specified in terms of the 4-vectors of the incom- 




Figure 1; Illustration of the generic process 7p 
exchanged between the photon and the proton. 



XY . In the Regge pole picture, a reggeon is 



ing photon and proton, denoted by q and P respectively, and those of the two distinct 
components of the final state, p^ and p^ . Convenient Lorentz scalars are 



W 



{q + Pf 



t={P-V, 



Mi 



v\ 



Mi 



vl. 



where W is the total photon-proton centre of mass energy, t is the square of the 
4-momentum transferred and M^ and My are the invariant masses of the two final state 
hadronic systems. The ratios M^/VT^ and IvP' jW'^ most naturally distinguish between 
the kinematic region where diffraction is dominant and that in which other exchanges 
become important. 

The mass of the system X is measured using the central components of the detector up 
to M^ ~ 30 GeV. For M^ > 10 GeV the mass of the system Y may also be measured. 
Where M^ is smaller its value is constrained using the location of the rapidity gap. The 
resolution in transverse momentum is insufficient for a differential measurement of the 
t dependence. Measurements of the differential cross section da/dM^ are presented at 
(PF) = 187 GeV and (PF) = 231 GeV, for the region M^ < 1.6 GeV and |t| < 1.0 GeV^ 
and for 1.6 < M^ < 15.0 GeV and all |t|. 



3 The HI detector 



The components of the HI detector that are most relevant to this analysis are briefly 
described here. More detailed information can be found in 



26[| . The coordinate system 



convention for the experiment defines the forward, positive z direction as being that of the 
proton beam. The polar angle 6 is defined relative to this axis such that pseudo-rapidity, 
1] = —In tan ^/2, is positive in the forward region. 

Final state electrons are detected in the electron tagger of the luminosity system (eTag), 
which is situated at z = —33 m and has an acceptance such that Q^ < 10~^ GeV^, where 
Q^ describes the photon virtuality. Bremsstrahlung photons produced at very small polar 
angles leave the beam-pipe through an exit window and are detected at 2; = —103 m in 
the photon arm of the luminosity system. 

Minimally biased photoproduction data samples are obtained using trigger signals from 
the eTag in coincidence with a 'ToF' trigger. The ToF trigger demands at least one 



6 



pair of hits in the twin layers of plastic scintillator comprising the Time of Flight device, 
during a short time period at each bunch crossing when particles originating from ep 
interactions in the vertex region are expected to arrive. The backward location of the 
ToF (—3.5 < 7] < —2.0) significantly enhances the efficiency for triggering events at low 
M^ compared to triggers that require tracks in the central parts of the detector. 

The system X is detected and measured using tracking and calorimeter information from 
the central parts of the detector. Charged tracks are measured in the range —1.5 < r] < 1.5 
in the two large concentric drift chambers of the central tracker, which lies inside a 
1.15 T solenoidal field. Surrounding the drift chambers, the finely segmented liquid argon 
calorimeter (LAr) provides smooth and hermetic coverage in the range —1.5 < rj < 3.4. 
The lead-scintillator Backward Electromagnetic Calorimeter (BEMC) completes the cov- 
erage in the backward region for —3.5 < f] < —1.5. The backward multi-wire proportional 
chambers (BPC) are mounted directly onto the surface of the BEMC. 

The most forward components of HI are used to tag hadronic activity at high pseudo- 
rapidity pTj-ppj]. At all but the smallest values of M^, rescattering of primary particles 
with the beam-pipe and the surrounding material yields showers of secondary particles 
which are observable in these detectors, extending their sensitivity to pseudo-rapidities 
larger than would be expected on the basis of their geometrical acceptance alone. The 
Proton Remnant Tagger (PRT) consists of double layers of scintillator surrounding the 
beam-pipe at z = 24 m and covers the highest pseudo-rapidity region. Three double layers 
of drift chambers of the Forward Muon Detector (FMD) cover a slightly lower pseudo- 
rapidity range. These components of HI provide high efficiency for the identification of 
primary particle production in the region 5 ^ rj < 7. A third detector component, the plug 
calorimeter, is sensitive at the lower end of this range and is used to perform cross-checks of 
the measurement. The overall range of sensitivity of the HI detector to final state hadron 
production spans 11 pseudo-rapidity units, providing excellent discrimination between 
different event topologies. 



4 Simulations 



Two Monte Carlo models, based on the event generators PHOJET [^ and PYTHIA |32 



are used to determine corrections to the data for triggering and detector inefficiency, the 
effects of analysis cuts and the smearing of kinematic quantities due to detector resolution. 
The description of soft photoproduction processes by both models is based on Regge 
phenomenology. Both generators produce events from all four elastic and dissociative 
subprocesses described in section ^ as well as events in which colour is exchanged. The 
generated events are passed through a full simulation of the HI detector and are subjected 
to the same analysis chain as the data. 

The lowest lying vector meson states are generated in the ratio p^ : u : (p = 14 : 1.5 : 1, 
consistent with previous measurements at HERA [p] -|23[| and at lower energy [Q. The 
differential t distributions take the peripheral form da^^^/dt oc e el*, with the slope param- 



eter in both models 6^^ — H GeV~^ for W = 200 GeV. Non- resonant di-pion production 
is accounted for by reweighting the mass distributions produced by the models such that 
the p° line shape matches the parameterisation of Ross and Stodolsky |3^ , with parame- 
ters fixed according to previous HERA measurements pT| , |2^ . 

Dissociative events are produced in the regions M"^ /W"^ < 0.1 and M"^ /W"^ < 0.1, using 
approximations to triple-pomeron predictions with couplings constrained by low energy 
data. The t-slope parameters, b^^, bp^ and 6^^, have logarithmic dependences on M^, 
My and W. The PD and GD slope parameters are approximately half those for the EL 
reaction, and the DD slope parameter is reduced by a further factor of approximately 2. 

The PYTHIA model of dissociation processes assumes a pomeron intercept of unity, such 
that the differential cross section for the GD process at fixed W takes the approximate 
form 

dtdM| M| ' ^ ^ 

with similar expressions for the differential PD and DD cross sections. Additional factors 
are applied in order to modify the distributions in kinematic regions in which a triple- 
pomeron model is known to be inappropriate. Their main effects are to enhance the 
low mass components of the dissociation spectra, suppress the production of very large 
masses and, in the DD case, to reduce the probability of the systems X and Y overlapping 
in rapidity space ||32|, ^]. GD, PD and DD states produced with M^ < 1.8 GeV or 
My < 1.8 GeV decay isotropically to two-body final states. The fragmentation of higher 
mass dissociation systems is treated in a string model, with final state hadrons distributed 
in a longitudinal phase space with limited transverse momentum p^ [^ . 



In the PHOJET model, the first moment Finite Mass Sum Rule (FMSR) ||3^ is applied 
in order to achieve a smooth transition between the elastic, resonance and high mass 
continuum regions. At fixed W, the GD spectrum follows the parameterisation 



r^ — 

dtdu yu 



— --(-) c^gdH* , (2) 



where the pomeron intercept is taken to be a]p(0) = 1.07, z/ = M^ — M^ — t and the 
'mass' Mq of the incoming photon is taken to be that of the p meson. A similar para- 
meterisation was found to work well for fixed target photoproduction data [PS. PD and 



DD distributions are produced in a similar manner with the proton mass defining M^ for 
the system Y. Dissociation to states in the resonance region (M^ < 5 GeV^) are treated 
through the extended Vector Dominance Model |3^ in a similar manner to EL processes. 



A single effective resonance is assumed in the low mass GD region, normalised using the 
FMSR and decaying to 2 or 3 particle final states. The slope parameter b,^^ falls smoothly 
through the resonance region. Higher mass final states are again limited in p^ . 

The mass distributions and ratios for the distinct subprocesses are iteratively reweighted 
such that the differential cross sections of each Monte Carlo simulation best describe the 



measured cross sections presented in section pTl . 



A third event generator, DIFFVM [BSI, models vector meson production at large Q^, 
both where the proton dissociates and where it remains intact. The simulation is based 
on Regge behaviour and the Vector Dominance Model and is described in more detail 
in [p9| ,pO[]. Since DIFFVM is able to generate dissociative events both with limited p^ 
and with isotropic decays of the dissociative system, it is used to check the sensitivity of 
the measurement to differing decay properties of the system Y. 



5 Experimental method 

5.1 Selection of photoproduction events 

The data used in this measurement were collected in two short dedicated runs in 1994 



when HERA was colliding positrons at 27.5 GeV with 820 GeV protons. The first run 
was taken with the vertex at its nominal position of z = 4 cm and the second with the 
vertex shifted in the proton direction to z = 71 cm. The latter configuration gives an 
enhanced efficiency for the detection of the most backward going final state particles. 
The integrated luminosities of the two samples are 24.7 ± 0.4 nb~^ and 23.8 ± 1.3 nb~^ 
respectively. The data samples have previously been used in measurements of the total 
7P cross section [EO] and of the cross section for exclusive p^ photoproduction 



Photoproduction events were triggered on the basis of a coincidence of signals from the 
ToF and eTag triggers described in section §. These requirements have been shown to 



give good efficiency for all photoproduction subprocesses ||20|. The eTag is also used to 
measure the energy of the scattered electron, from which the 7p centre of mass energy W 
is inferred. Two ranges of W are considered, covering the region in which the acceptance 
of the eTag is highest. The first, 164 < W < 212 GeV, yields a centre of mass energy 
averaged over the photon fiux distribution 0,^ of {W) = 187 GeV. The second range 



is 212 < VT < 251 GeV, for which (W) = 231 GeV. 

Events arising from the Bremsstrahlung process, ep -^ eprf, are removed by requiring no 
activity above 2 GeV in the photon detector. A correction of 1.8% is made for the loss of 
signal due to the random overlap of Bremsstrahlung events with '-/p interactions. 

The principal background to the measurement arises from interactions of the electron 
beam with residual gas particles or with the walls of the beam-pipe. Proton-beam induced 
interactions also contribute to the background when they occur in coincidence with a 
reconstructed electron in the eTag. This may arise from electron-gas/wall collisions or 
from a Bremsstrahlung process in which the photon is not detected. All of these sources 
of contamination are suppressed by requiring either that the total energy reconstructed in 
the calorimeters is greater than 1.5 GeV, or that a vertex is reconstructed at a position in z 
that lies within 30 cm (~ 3cr) of the mean interaction point. Any event in which a vertex is 
reconstructed outside this region is rejected. In the higher range of W, where backgrounds 
are largest, a BPC hit is required in addition. After these cuts, the remaining level of 



background associated with the electron beam alone is estimated from non-interacting 
electron pilot bunches and that arising from the proton beam is determined from event 
samples in which either the ToF or the eTag component of the trigger is relaxed. Both 
backgrounds are subtracted statistically. The subtraction varies between zero and 19% in 
the measured region. 



5.2 Separation of the systems X and 1^ 

In order to ensure that the systems X and Y are clearly separated, only those events 
that contain a region of at least two units of pseudo-rapidity in which no particles are 
reconstructed are considered in the analysis. This gap extends from the largest pseudo- 
rapidity reached by fragments of the system X to the smallest pseudo-rapidity at which 
particles associated with the system Y are observed. Its presence implies that both M^ 
and My are small compared to W. 

Three distinct subsamples are defined. In sample A there is no activity above noise levels 

in the forward detectors (PRT and FMD) and the most forward part of the LAr {6 < 4.7°). 

These cuts are similar to those which have been used in previous HI publications 

2^, p7|-|30|. With this selection, the pseudo-rapidity gap spans at least 3.2 ^ rj < 7.0 



under nominal vertex conditions and 3.0 ^ rj < 7.0 for shifted vertex data. In sample 
B a signal is observed above noise levels in the forward detectors and there is a pseudo- 
rapidity gap of at least two units in the central detectors (central tracker and LAr and 
BEMC calorimeters), extending to the forward limit of acceptance of the LAr. Sample C 
is defined in the same way as sample B, except that the gap in the central detectors does 
not extend to the forward edge of the LAr. 

The data and simulated events that are shown in all subsequent figures are required to 
pass one of these three selection criteria. There are 8685 such events in the nominal vertex 
data sample and 8184 in the shifted vertex sample. 



5.3 Reconstruction of M^ and M^. 

The mass reconstruction has been optimised by combining tracking and calorimeter in- 
formation in an energy flow algorithm that avoids double counting. Primary vertex con- 
strained tracks are extrapolated into the calorimeter and energy clusters increasingly 
distant from the extrapolated track are discounted. The procedure continues either until 
the total energy of the excluded clusters exceeds that of the track, under the assumption 
that the particle yielding the track is a pion, or until all electromagnetic (hadronic) clus- 
ters within a cylinder of radius 30 cm (50 cm) about the extrapolated track have been 
removed. Detector noise is potentially problematic, particularly at small values of M^^. To 
minimise its effect, techniques have been developed that reject low energy isolated clusters 
in the LAr. Any calorimeter deposit that is not associated with a track or rejected as 
noise is accepted. 
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The same method is used to reconstruct M^ for all three samples A, B and C The gap 



requirements introduced in section ^]2| ensure that X does not extend beyond the forward 
acceptance limit of the central detector components. Losses in the backward direction 
are, however, unavoidable both into the beam-pipe and because of the poor containment 
of hadrons by the BEMC. In order to minimise the effect of such losses, the mass of the 
system X is reconstructed using the expression 



M' 



2 q.p^ +Q^ + t ~ 2E^ Y.^E + p. 



(3) 



where E^ is the laboratory energy of the interacting photon, obtained from the electron 
energy measurement, and t and Q^ have been neglected in the last expression. The 
quantity Y.{.E +Pz)x i^ summed over all hadrons reconstructed backward of the largest 
pseudo-rapidity gap and has the property of being insensitive to the very backward going 
hadrons with E ~ —pz that are usually not observed. A constant correction factor of 1.07 
is applied to the reconstructed M^ to account for the remaining losses. Figure ^ shows 
a comparison between the reconstructed and generated values of M^ in the PYTHIA 
simulation. Good reconstruction is obtained for all masses, including those corresponding 
to exclusive vector meson production. The resolution shows a slow improvement with 
increasing mass and is approximately 30%. 
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Figure 2: Correlation between reconstructed and generated values of (a) M^ and (b) M^ for events 
in the shifted vertex PYTHIA simulation with 164, < W < 251 GeV. 

Only in sample C, where fragments from the system Y are observed in the central detector 
components, is it possible to reconstruct M^ ■ For this sample a similar approach is taken 
for the reconstruction of M^ as for M^ . The relation 



M' 



2p.p^ 



m 



t ~ 2Ep Y.iE 



P. 



ZJY 



(4) 



is used, where E^ is the energy of the proton beam and the proton mass m and t are 
neglected in the last expression. Equation (H) is insensitive to losses at very large pseudo- 
rapidities. The remaining effect of undetected particles is accounted for by applying a 
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constant correction factor of 1.10 to the reconstructed My. Figure 0b shows a comparison 
between reconstructed and generated My values from simulated events in which particles 
from the system Y reach the central detector components (sample C). Good reconstruction 
is obtained for M, > 10 GeV. 



Y 



For samples A and B the system Y is not observed in the central detectors and direct 
reconstruction of My is not possible. Instead, information is obtained from the presence 
or absence of activity in the forward detectors. The location of the upper edge of the 
pseudo-rapidity gap is correlated with My so that, with a knowledge of the regions of 
pseudo-rapidity to which the forward detectors are sensitive, different regions of My are 
distinguished. 



5.4 Selection of intervals in the invariant mass M^ 



The three subsamples defined in section ^^ are used to measure cross sections in two 
ranges of M^. For sample A the system Y passes unobserved down the forward beam 
pipe at r] > 7. Studies using the Monte Carlo simulations have shown that this selection is 
almost 100% efficient for events in which y is a single proton and falls to approximately 
50% at My = 1.6 GeV. Sample A is therefore used to measure the cross section with 
My < 1.6 GeV. This region is likely to be dominated by the subprocesses GD and EL 
in which y is a proton. There are also likely to be contributions from events in which Y 
is a different low mass baryonic system. Since neither the isospin nor the charge of Y is 
determined, this may be any of n, N*, A, . . ., or a non-resonant system. 

Samples B and C are used together to measure the cross section for events of the types DD 
and PD integrated over the range 1.6 < My < 15.0 GeV. For sample B, My < 5 GeV, 
while in sample C My is usually larger. The full sample B is used along with events 
from sample C for which My, reconstructed using equation (|), is less than 15.0 GeV. 
The restrictions imposed by the pseudo-rapidity gap criteria mean that the maximum 
accessible value of M^ decreases as My increases. The upper limit of My = 15.0 GeV is 
chosen so that the measurement can be made over a large range of M^ , whilst ensuring 
that the upper limit in My lies in the region where the reconstruction is good. 

Figures ^ and ^ show uncorrected M^ distributions for events used to make measure- 
ments in the two My intervals compared to the simulations. Satisfactory agreement 
between the data and the Monte Carlo models is obtained in both cases. An elastic 
peak is observed in sample A, along with a higher mass continuum that extends until the 
system X reaches the limit set by demanding no activity in the most forward region of 
the LAr. In samples B and C no significant enhancement in the logM^ distribution is 
observed in the region of the light vector mesons, and the requirement of a gap of at least 
two pseudo-rapidity units in the central detector components restricts the data to smaller 
values of M^ . Figure ^ shows the reconstructed distribution in proton dissociation mass 
of events from data sample C. 
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Figure 3: Uncorrected mass distributions for the shifted vertex data sample in the range 
164 < W < 212 GeV compared to the two simulations normalised via the total 7p cross section |g^. (a) 
M^ distribution of events from sample A, used to measure the cross section for the region M^ < 1.6 GeV. 
(b) As in (a), but for events in samples B and C that are used to measure the cross section for 
1.6 < M^ < 15.0 GeV. (c) Distribution in M^ of events in sample C, where activity from Y is ob- 
served in the central detector components. 

5.5 Extraction of cross sections 



The criteria described in section ^.4| are used to measure cross sections as defined in 
section H for tfie regions M^ < 1.6 GeV and 1.6 < My < 15.0 GeV. Tlie events in eacli M^ 
range are divided into intervals in M^ and tlie residual contribution from beam induced 



background is subtracted in each interval using the procedure described in section 5.1 



An average of the PHOJET and PYTHIA models is used to correct for losses and mi- 
gration of events in M^ and My. For |t| > 1.0 GeV^, elastically scattered protons are 
observed in the PRT, resulting in their classification in the wrong M^ interval. The sim- 
ulations are therefore also used to correct the measurement for My < 1.6 GeV to the 
region |t| < 1.0 GeV^ and that for 1.6 < M^ < 15.0 GeV to all |t|. 

The Monte Carlo models do not explicitly simulate the states Y that may be present in 
the data but have different isospin or charge from the proton. However, the acceptance for 
such states is likely to be similar to that for final states arising from diffractive exchange 
and, after the reweighting procedure described in section |^, no additional corrections are 
made. A small correction is applied for the loss of events contributing to the cross section 
at large My and large M^ for which the largest gap in pseudo-rapidity is smaller than 
two units. The correction necessary in unfolding from the experimental separation of the 
systems X and Y using the largest gap in pseudo-rapidity to the cross section defined in 
section ^ in terms of rapidity is found to be negligible. 

The mean acceptance of the eTag is 64 ± 2% for the (W) = 187 GeV sample and 35 ± 2% 
for (W) = 231 GeV. The combined acceptance of the ToF component of the trigger 
and the off-line cuts is greater than 40% for all measurements and purities^] are better 

^Purity is defined as the fraction of all simulated events reconstructed in a mass interval that are also 
generated in that interval. 
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than 30% for all measured data points. The maximum values of M'^/W'^ and M'^/W'^ 
considered are 0.041 and 0.0084 respectively. Migration from the regions M"^ /W'^ > 0.1 
and M"^ /W"^ > 0.1 is less than 10% for all quoted measurements. 

A 7p cross section is obtained assuming a photon flux derived from the Weizsacker- 
Williams approximation pO| , ^ . Radiative corrections were found to be at the level of 1% 



in [^ and are neglected here. The measurements obtained from the shifted and nominal 
vertex data samples are compatible within statistical errors. The two measurements are 
averaged with weights determined from their statistical errors to produce the quoted cross 
sections. All measurements are presented at the bin centres in log M^ . With the exception 
of the lowest two M^ intervals, where exclusive vector meson production gives rise to a 
resonance structure p2|, bin centre corrections at the level of 2% are applied, based on 



the fits presented in sections q]3| and pA . 



For (W) = 231 GeV the boost of the jp centre of mass frame relative to the laboratory is 
reduced compared to that at (W) = 187 GeV, such that measurements in the lowest M^ 
region are no longer possible, but larger M^ values become accessible. Statistical errors 
are largest at high W, since both the photon flux distribution and the acceptance of the 
eTag fall with increasing photon energy. 

A cross check has been performed with a sample of data triggered by a track in the 
central chambers [Q, for which the acceptance becomes high at large M^. The results 



were consistent with those obtained from the ToF triggered samples. The sensitivity of the 
results to variations in the pseudo-rapidity gap selection criteria and mass reconstruction 
methods have been checked in several ways. No significant changes in the measured cross 
sections are observed when the forward region of the LAr that must be devoid of activity in 
the low My measurement is enlarged to 6' < 10° or when this cut is removed completely. 
There is also little change when the minimum gap size requirement for the larger My 
measurement is varied between 1 and 3 pseudo-rapidity units, when a further forward 
detector, the plug calorimeter, is included, or when M^ and My are reconstructed using 
calorimeter clusters only. 



5.6 Error analysis 

Statistical errors arise from the finite volume of data, the samples used in the subtraction 
of beam induced background and the simulation samples used to make corrections. The 
statistical error in each interval is formed by adding the contributions from these three 
sources in quadrature. 

Systematic errors are evaluated on a bin by bin basis. No single contribution dominates 
in any region of the measurement. A brief summary of the systematic effects that are 
found to be most important is given below. 

• The uncertainties in the hadronic energy scales of the BEMC and LAr calorimeters 
are 20% and 5% respectively. That on track measurements is 3%. There is a 1% 
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uncertainty in the energy calibration of the eTag. 

• The uncertainty in the efficiency of the PRT is 20%. 

• A 10% variation in the beam induced background, accounting for uncertainties in 
the extrapolation of pilot bunch information. 

• Uncertainties in the decomposition of the total photoproduction cross section are 
accounted for with 50% variations in the subprocess ratios (y^^/cjJot and Cp^/o^Zot in 
the simulations and 20% variations in the ratios a^^/a'^^^^ and a^^/alot- 

• The shapes of the GD and DD mass spectra assumed in the simulations are indi- 
vidually reweighted by a factor {^r)^'^'^^- Those of the PD and DD processes are 

reweighted by a factor {^r)^^'^- 

• The simulated t dependence is varied by changing the slope parameters by ±4 GeV~^ 
for EL, ±2 GeV-2 for GD and PD and ±1 GeV'^ for DD processes. 

• A 50% variation is applied in the relative strengths of uo and production in the 
simulations. 

• Uncertainties in the correction for migration between the two My intervals are 
estimated by varying the Y fragmentation scheme assumed in the simulations. An 
error is formed from the spread in corrections obtained using PHOJET, PYTHIA 
and DIFFVM in isotropic and longitudinal phase space modes. 

• The full difference between results obtained when correcting the data with the PHO- 
JET and PYTHIA simulations is taken as the uncertainty in the modelling of the 
final state. 

• The errors on the luminosity measurements result in overall normalisation uncer- 
tainties of 1.7% for the nominal vertex samples and 5.4% for the shifted vertex 
samples. 

• The uncertainty in the acceptance of the eTag results in overall normalisation un- 
certainties of 3.5% averaged over the range 164 < W < 212 GeV and 6.0% for 
212 <W < 251 GeV. 

The systematic errors quoted for each data point are formed by adding in quadrature the 
variations in the measured cross section arising from each uncertainty. 



6 Results 
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6.1 The differential cross sections dcr/dM 



The measured differential cross sections are presented in figure ^. Figure ^ shows 
dcr^p^xy/dM^ at (W) = 187 GeV, measured in 9 bins of M^ in the range 
0.160 <Ml< 862 GeV^ with M^ < 1.6 GeV and |t| < 1.0 GeV^. Superimposed is 
the cross section integrated over the region 1.6 < M^ < 15.0 GeV and all |t|, in 6 M^ 
intervals in the range 0.160 < M^ < 86.2 GeV^. The first interval in M^ is chosen to 
cover the region of p", uj and (p production. The second is constructed to contain the p' 
For these two measurements an error bar, reflecting the width of the M^ 



resonances 



41 



interval, is shown to account for uncertainties in the details of the resonance structure.^ 
For all but the first two M^ intervals measurements are made in three bins per decade 
of M^ . The results of the measurement at (W) = 231 GeV are shown in figure ^. All 
cross section measurements are summarised in table |I[ 
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Figure 4: The measured differential cross section dcr/dM^ for the process 7p -^ XY in the regions 
M^ < 1.6 GeV, \t\ < 1.0 GeV2 and 1.6 < M^ < 15.0 GeV. (a) (TV) = 187 GeV, (b) (W) = 231 GeV. 
The inner error bars show statistical uncertainties only. The outer ones show statistical and systematic 
errors added in quadrature. 



At sufficiently large M^ the measurements for My < 1.6 GeV in both W ranges dis- 
play an approximate dcr/dM^ ~ 1/^^ dependence, though closer scrutiny reveals a 



more complex structure which is investigated in section |6.3| . The M^ distribution for 



•^In the remaining M^ intervals the largest deviation from a smooth dependence on M^ is likely to 
arise from J/^ production, which contribute at the level of 5% in the fourth M^ interval |3y] . 
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(a) 







{W) = 187 GeV 


{W) ^ 231 GcV 


Ml 
(GeV^) 


Bin Limits 
(GeV2) 


g|^ Stat. Syst. 
"" {^ih GeV-2) 


^ Stat. Syst. 
" [ph GeV-2) 


0.504 


0.16-1.58 


(1.25 ±0.05 ±0.11) xlO^ 


— 


2.52 


1.58-4.00 


(1.80 ±0.11 ±0.22) xlO" 


(1.75 ±0.26 ±0.24) xlO*^ 


5.87 


4.00 - 8.62 


(5.37 ±0.34 ±0.48) xlO"^ 


(6.67 ±1.00 ±0.92) xlO-^ 


12.7 


8.62-18.6 


(2.24 ±0.15 ±0.22) xlO"^ 


(2.13 ±0.32 ±0.24) xlO"^ 


27.3 


18.6-40.0 


(8.25 ±0.74 ±0.52) xlO"^ 


(8.68 ±1.76 ±1.12) xlO"^ 


58.7 


40.0-86.2 


(3.53 ±0.34 ±0.25) xlO"^ 


(4.07 ±0.72 ±0.36) xlO^^ 


127 


86.2-186 


(1.79 ±0.15 ±0.19) xlO-^ 


(2.10 ±0.28 ±0.21) xlO-^ 


273 


186 - 400 


(1.01 ±0.76 ±0.91) xlO"^ 


(1.05 ±0.14 ±0.10) xlO""^ 


587 


400 - 862 


(5.56 ±0.40 ±0.55) xlO"^ 


(5.50 ±0.65 ±0.73) xlO"^ 


1270 


862 - 1860 


— 


(3.22 ±0.48 ±0.52) xlO"^ 



(b) 









w = 


187 GeV 


{W) = 231 GcV 


Ml 
(GeV^) 


Bin Limits 
(GeV^) 




Stat. 
[^lh GeV 


Syst. 


da 


Stat. Syst. 
;^b GeV-2) 


0.504 


0.16- 1.58 


(3.40 


±0.25 


±0.40) xlO'J 


— 


2.52 


1.58-4.00 


(7.60 


±0.80 


±1.38) xlO-^ 


(8.91 


±2.34 ±1.34) xlO"^ 


5.87 


4.00-8.62 


(2.98 


±0.34 


±0.49) xlO-i 


(2.59 


±0.79 ±0.46) xlO-i 


12.7 


8.62-18.6 


(1.79 


±0.18 


±0.30) xlO-^ 


(1.33 


±0.35 ±0.20) xlO-^ 


27.3 


18.6-40.0 


(8.19 


±0.82 


±1.44) xlO"^ 


(8.57 


±1.79 ±1.22) xlO"^ 


58.7 


40.0-86.2 


(4.56 


±0.51 


±0.85) xlO"^ 


(3.51 


±0.84 ±0.79) xlO"^ 



Table 1: The measured differential cross section dcr/dAf^ for the ranges (a) M^ < 1.6 GeV and 
|i| < 1.0 GeV^ and (b) 1.6 < M^ < 15.0 GeV and all t. All measurements are quoted at the bin centres 
in logM^ . Overall scale uncertainties of 5.2% at (W) = 187 GeV and 6.9% at (W) = 231 GeV are not 
included in the systematic errors. 



^ than that for 



1.6 < My < 15.0 GeV falls slightly less steeply with increasing M 

My < 1.6 GeV. The cross section is enhanced considerably at the lowest values of 

Ml and M^ where elastic processes dominate 12 

Previous data exist on the differential cross section for 7p -^ Xp ||25]. Since the region 



My < 1.6 GeV is expected to be heavily dominated by this process, the low M^ mea- 



surement is compared to fixed target data in section |6.3| and combined fits are performed 
using the Regge pole model described in section |0| . The degree to which diffractive and 
non-diffractive exchanges contribute to the single photon dissociation process is investi- 
gated. The questions are addressed as to whether the intercept of the pomeron in such 
a model is consistent with that previously deduced from soft hadron-hadron and pho- 
toproduction interactions and whether any anomalous behaviour is necessary to explain 
the W dependence of the diffractive contribution. The behaviour of the cross section for 



1.6 < My < 15.0 GeV is discussed in section p^ . 
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6.2 Triple-Regge model of photon dissociation 

An appropriate framework in which to model dissociation processes is offered by 'triple- 
Regge' phenomenology. For W"^ ^ M^ the process 7p -^ Xp may be treated with a Regge 
expansion. The amplitude at fixed M^ is then a sum of amplitudes for the exchange of 
reggeons i that produce all possible final states X, as illustrated in figure ^. The corre- 
sponding cross section contains products of flux factors for reggeons i and j and amplitudes 
T^ai{t)^x ^^Q,.(j)^x (figure ^d). The generalised optical theorem [0 relates the sum over 
X of these matrix elements to the forward amplitude for the process 70:4 (t) -^ lotji^) 
at an effective centre of mass energy M^. When M^ is large by comparison with the 
hadronic mass scale Sq, a Regge expansion is also appropriate for the photon-reggeon 
scattering amplitude, such that the dissociation cross section may be decomposed into 
triple-Regge terms as shown in figure ^. The cross section may then be expressed as a 
sum of contributions with reggeons z, j and k [§-[11]] 
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Figure 5: Illustration of the Mueller- Regge approach to the inclusive photon dissociation cross section. 

The functions Gijk{t) and ai{t) are not predicted by the model and must be determined 
from experimental measurements. The trajectories aj(t) are assumed to take the linear 
form aiit) = ai{0)+a'^t. The phase 0j(t) of reggeon i is determined by the signature factor, 
r]i{t) = (^ -\- e~™°'^^^\ where C. = ±1 is the signature of the exchange. The signature factors 
are written as i]i(t) = ril{t) e*"'^''-*^ with the moduli rfl{t) absorbed into the (3 parameters 
introduced in equation (^. For photoproduction reggeons i and j must have the same 
signature such that 0j(t) — (pjif) = | [aj{t) — ai(t)]. As is customary the scale s^ is set 
to 1 GeV^. The functions Gijkit) may be factored into products of couplings each of the 
forniQ 



Gijkit) 



167r 



f^pi{t) Ppjit) P^kiO) Qijkit) 



(6) 



"'The couplings Gijk{t) are expected to differ from those extracted from pp data by a factor 

(3^ki0)/(3pki0). 
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where the j3 terms describe the couphngs of the reggeons to external particles and 
gijkif) is the appropriate three-reggeon coupling. The t dependence of the reggeon- 
proton and three-reggeon couplings are parameterised here as (3pi{t) = /5pj(0)e''p'* and 
9ijk{t) = 5'ijfc(0)e^>^'=*. 

The pomeron, with trajectory dpit), is unique in Regge theory in that its intercept is 
significantly larger than those of all other reggeons. In the limit in which both W'^/M'^ 
and M^ -^ oo only the ijk = PPP term survives and equation (Bl) reduces to 



X 



da G'ppp(O) ^^2\2«p(0)-2 /J_y^^^^ ^B{W^Ml)t /7^ 



X oo \ X / 

where B{W'^, M^) = 26pp + 6ppp + 2a'^ In W'^/M'^. After the pomeron the next-leading 
reggeons have approximately degenerate trajectories and carry the quantum numbers of 
the p, uj, a2 and /2 mesons.0 In this analysis the symbol IR is used to describe combinations 
of these four reggeons and a single effective trajectory a^^it) is assumed. 

With the two trajectories, apit) and ^^^(t), equation @ leads to a total of six terms 
with distinct W^ and M^ dependences. Table ^ lists each of these terms. Diffractive 
contributions are approximately independent of W"^ and correspond to the case in figure |^c 
where both reggeons i and j are the pomeron. In addition to the triple-pomeron diagram, 
a further diffractive term arises from the ijk = PPM diagram. The reggeon k must have 
the quantum numbers of the / in order to satisfy the requirements of conservation of 
C-parity at the photon vertex and C- and G-parity at the three-reggeon vertex. 

The two non- diffractive terms, PPP and PPP, may involve the exchange of any of p, 
u, a or f and are suppressed by factors like l/W"^ at fixed M^ . The production of states 
Y that have different isospin or charge from the proton {Y = n, N*^, A"'"^, A"*", A°, . . .) 
can only occur via isovector exchanges {i = j = p or a). Such processes are described, 
up to a normalisation, in the same way as processes that are elastic at the proton vertex, 
through those specific terms in table ^ marked with a star. 

Since the P and / reggeons exchange identical quantum numbers they may inter- 
fere [^^ -^.pl Such a contribution is expressed in terms of the 'non-diagonal' amplitudes, 
{PP}P and {PP}1R where {PP}fc = JPJRk + JRJPk. The / reggeon has previously 



been found to couple strongly to the proton [^ and since the interference terms have 
weaker dependences on W'^ than the non-diffractive contributions (PPP and PPP), 
they may be of particular significance in the high energy regime of the present measure- 
ment. Restrictions are imposed on the magnitudes of the interference relative to the 



''The reggeons under consideration are hereafter referred to as P, p, u, a and /. Their isospin, 
signature and C- and G-parities are P(0+++), p(l h), w(0 ), a(l +H — ) and /(0+++). 

^Only interference between the W and / reggeons is considered in this analysis. Since the photon is not 
an eigenstate of G-parity, interference between reggeons of opposite G-parity is in principle also possible. 
However, the dominance of two-pion final states at small M^ suggests that the photon is dominated by 
G = -t-1 and that the interference terms JP aa, poja and /oo are likely to be small [EtI. 
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non-interference terms through the inequahty ^ 



=^"*'°'EG...w(Mr"''°' 



Generic 
term ijk 


Specific 
terms ijk 


Approx. W^ 
dependence 


Approx. M^ 
dependence 


PPF 


PPP 


lyu 


M- 


FPU 


PP/ 


VK" 


M-'-^ 


PIRP 


//P aaP* ppP* wwP 


VK-^ 


Ml 


PPP 


/// aa/* PPf* ^^f {/a}a {p^ja 


ly-^ 


M- 


{PP}P 


{P/}P 


w^-1 


M;i 


{PP}P 


{P/}/ {Pa} a 


w~' 


M- 



Table 2: Tire VK and M^ dependence of eacli of the six triple- Regge amplitudes considered in the fits 
under the approximations a-p{t) ~ 1 and a-^{t) ~ 0.5. The second column lists each of the particular 
terms that are allowed by quantum number conservation. Those marked with a star can result in states 
Y that have different isospin or charge from the proton. 



The intercept of the pomeron is a free parameter in the fits presented here. A free fit 
for all other parameters is not practical with the limited photoproduction data available. 
Instead, information is taken from data on total and elastic cross sections. The intercept 



of the subleading trajectory is taken from the results of fits |24[ to the centre of mass 
energy (y^) dependence of total cross sections that use expressions of the form 



^ab _ j^ab^a^(Q)-l _|_ Qab ^a^{0)^ 



(9) 



The trajectory gradient a' is obtained from the shrinkage of the forward peak in high 
energy elastic pp scattering [|^. The slopes for all non-diffractive trajectories are assumed 



to be the same and the value of a'^ is taken from the result obtained for a' in analyses [^ 
of the reaction 7r~p — > vr^n. The slope parameter hp^ is determined from elastic scattering 



measurements [49| in the regime of diffractive dominance using 



dt 






M^)i 



(10) 



t=0 



where B{s) = 46pp + 2aMns/so. The parameter bp^^i is less well determined and is 
allowed to vary within a range that covers most of the results obtained in lower energy 



fits. Previous measurements [£T| have found that the t-dependence of the triple-pomeron 
coupling is weak compared to those of the external vertices, though less is known about the 
remaining three-reggeon couplings. It is assumed here that the dominant t dependence of 
Gijk{t) arises from the (3 terms and uncertainties in the t dependence of the three-reggeon 
couplings are accounted for by setting bijk = ± 1 GeV~^ for all ijk. Table ^ specifies 
the values assumed for the parameters that are kept fixed in the fits and the assigned 
uncertainties. 
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Quantity 




Value 


Source 


«m(0) 


0.55 


±0.10 




24 




< 


0.26 


±0.02 GeV-^ 




49 




< 


0.90 


±0.10 GeV-^ 




50 




bpjp 


2.3 


±0.3 GeV-^ 




49 




bpTR 


1.0 


±1.0 GcV-^ 


n 


bijk 


0.0 


±1.0 GeV-^ 





Table 3: Values assumed for the fixed parameters in the triple- Regge fits and the sources from which 
they are taken. Each quantity is defined in the text. Contributions to the model dependence errors are 
formed by repeating the fits after separately varying each parameter by the quoted uncertainties. 

6.3 Regge analysis of the proton-elastic cross section 



Large ranges in both W"^ and M^ are required in order to investigate the importance of all 
six terms in table |^. To extend the range in W"^ simultaneous fits are performed to HI and 
to fixed target data |^, where the quantity ^JdM^ ^^^ presented at t = —0.05 GeV^ 
in the ranges 11.9 <W< 13.7 GeV and 13.7 <W< 16.7 GeV. For the data from ^ 
equation (j^) is used at t = —0.05 GeV^. For the HI data equation (||) is integrated over 
the measured range, |tjnin| < N < 1 GeV^, where |tjnin| is the minimum kinematically 
accessible value of |t| |^. The parameters listed in table || are fixed and 0^(0) and 
the four couplings, G]pjpjp(0), G']p]p]j^(0), Gj^]j^j:,(0) and GnRMilO)) ^^^ ^^e fit parameters. In 
order to avoid the complications of the resonance region only data for which M 
are included in the fits. 



2 > 4 GeV2 



The additional contribution in the HI data from events in which the proton dissociates 
with My < 1.6 GeV is accounted for by applying a factor of 1.10 ±0.06 to equation (H) in 
the triple-Regge analysis for the HI measurements onlyQ. This correction is determined 
by taking an average of the values expected in the PHOJET and PYTHIA models after 
weights have been applied which make the differential cross sections in the models agree 
with those measured (see section |^). The uncertainty in this correction is taken as the 
difference in the values obtained from the two models. The correction does not account for 
additional states Y that may be produced by isovector exchange, but cannot be produced 
by diffractive or other isoscalar exchange. Any contribution from these states would 
increase the couplings Gj^jj^p(O) and G'jj^j^j^(O) for the HI measurement relative to the 
fixed target data (see table ||). 

No good solutions are found when diffractive terms only are considered. In a fit in which 
the pomeron alone contributes, equivalent to equation (|^) with a^p^O) and G^pp as free 
parameters, the x^/ndf made from statistical errors is 89.9/33. When both terms PPP 
and PPIR are included, the result is x^/^df = 79.3/32. We conclude that significant 
non- diffractive contributions are needed to explain the measured cross sections. 

^No correction is necessary for the fixed target measurements, which were performed by tagging 
leading final state protons. There may be a contribution from protons arising from decays of A and other 
resonances, though this effect is expected to be negligible p5|]. 
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(a) xVndf = 28.8/30 













Correlation 


coefficients 




Parameter 


Value 


Stat. 


Syst. 


Model 


^^ IP IP IP 


^^IPIPIR. 


^^IRIRIP 






a^(0) 


1.031 


±0.014 


±0.013 


±0.009 


-0.93 


-0.86 


-0.74 


0.82 




G'n.n.n, (0) 


8.19 


±1.60 


±1.34 


±2.22 ^bGeV-^ 




0.63 


0.50 


-0.60 




Gipipm^O) 


14.0 


±5.2 


±8.1 


±2.2 ^bGeV-^ 






0.89 


-0.94 




G„n..(0) 


238 


±73 


±98 


±101 ^bGeV-^ 








-0.98 




Gm,iRiR(0) 


-506 


±293 


±325 


±328 ^bGeV-" 
















(b) xVndf 


= 19.9/30 










Correlation coefficients 




Parameter 


Value 


Stat. 


Syst. 


Model 


IP IP IP 


^^IPIPIR. 








a^{0) 


1.101 


±0.010 


±0.022 


±0.022 


-0.88 


-0.46 


0.23 


0.08 




Gipipip(O) 


2.05 


±0.44 


±1.26 


±1.21 ^bGeV-' 




0.06 


-0.59 


0.32 




Gs>wm.{0) 


4.19 


±0.76 


±1.36 


±1.70 AibGeV"^ 






0.42 


-0.64 




CjRIRIP (0) 


115 


±29 


±65 


±78 ^b GeV"^ 








-0.93 




G'm,iR]R(0) 


-405 


±217 


±382 


±464 ^bGeV-" 
















(c) xVndf 


= 18.4/29 










Correlation coefficients 


Parameter 


Value 


Stat. 


Syst. 


Model 


IP IP IP 




^-^IRIRIP 




n 


a^{0) 


1.071 


±0.024 


±0.021 


±0.018 


-0.98 


-0.93 


-0.85 


0.92 


0.80 


Gipipip(O) 


3.76 


±1.62 


±2.25 


±1.45 AibGeV-^ 




0.83 


0.83 


-0.89 


-0.82 


Gipipin.(0) 


7.46 


±3.73 


±2.81 


±2.67AibGeV-^ 






0.84 


-0.90 


-0.75 


Cmmip (0) 


63.8 


±70.7 


±17.7 


±39.7AibGeV-^ 








-0.99 


-0.98 


Gm.iRiR(0) 


264 


±344 


±27 


±159 ^bGeV^'^ 










0.95 


7^ 


4.56 


±3.60 


±1.74 


±1.41 













Table 4: Results ol the triple-Regge fits and statistical, systematic and model related errors. The 
X^/ndf values and correlation coefficients reflecting statistical errors only are also given, (a) Fit with- 
out interference or isovector exchange contributions, (b) Maximal constructive interference between the 
diffractive and the secondary exchange, (c) No interference, but with a possible contribution from isovec- 
tor exchange. In fit (c) the sum of the couplings Gj^iRn=(0)+G'jR^j^j^(0) is allowed to be larger in the HI 
data than in the fixed target data by the factor TZ. The couplings shown pertain to the fixed target data. 



Three scenarios are examined in which subleading as well as diffractive contributions are 
considered. In fit (a) it is assumed that there is no interference between the diffractive 
and non- diffractive exchanges and that there is no contribution from isovector reggeons. 
Only the first four generic terms in table ^ are included and the IRIRP and HIRIR terms 
have the same normalisation at all W"^ values. Such a situation would be expected to 
occur if the uj were the dominant subleading reggeon. The fit parameters are shown in 
table |a. All experimental uncertainties are combined to determine the quoted systematic 
errors. The model dependence errors arise from the uncertainties in the parameters listed 
in table ^ and from the uncertainty in the correction from My < 1.6 GeV. All systematic 
and model dependence uncertainties are estimated from the differences between the main 
fit results and those obtained with the appropriate shifts imposed. 

Interference is introduced in fit (b) by including the {P1R}P and {PlRjlR terms. Fol- 
two-component duality is assumed in the reggeon-photon scattering 



lowing |45 
amplitude 



6,53 



such that the diagonal and non-diagonal terms are separately related for 
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Figure 6: Measurements of the quantity A'P da-yp^xy/dlVP with M^ < 1.6 GeV and |t| < 1.0 GeV^ 
by HI and of M^ da-yp^xp/dM^ di from |2^. For the HI data the inner error bars are statistical and the 
outer error bars show statistical and systematic errors added in quadrature. Overall scale uncertainties 
of 13% at (W) = 12.9 GeV and (W) = 15.3 GeV, 5.2% at (W) ^ 187 GeV and 6.9% at (W) = 231 GeV 
are omitted from the errors. The triple- Regge fit (b) with maximal constructive interference and the 
resulting decomposition of the cross section is superimposed. 



k = F and k = ]R 
interference couplings at t 
and G^^^j^, 



are parameterised as G. 



Maximal coherence is assumed and, through equation 

{IPIR.}]P V^/ 



the 



2\/|G'n,]p]p(0) Gj 



(0) = 2J|C^(0) G, 



This represents the scenario in which the / is 



the dominant subleading exchange and the / and the pomeron couple similarly to the 
proton [^Q. The results are presented in table ||b and figure |. 



In fit (c) effects arising from isovector exchanges are investigated. The interference terms 
are not included, but the non-diffractive terms (IRIRIR and IRHP) are allowed to be 
different in the HI and the fixed target data. This accounts for possible additional con- 
tributions to the HI data from the specific terms in table ^ that are marked with a star. 
A further free fit parameter TZ is therefore introduced, defined as the ratio of the sum of 
the couplings G'^]fj]p(0)+G'j^]j,]pj (0) in the HI data to that in the fixed target data. If there 
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were no isovector contribution to the subleading exchanges, then it would be expected 
that 71=1. The presence of p or a exchanges would increase 7^. The results of the fit 
are presented in table Ijc. 

From table ^ it can be seen that all three fits are acceptable. This demonstrates that the 
triple Regge model gives a satisfactory description of both the W^ and the M^ dependence 
of the existing data on photon dissociation, with the pomeron intercept in excess of 
unity. The x^/ndf value decreases markedly in fits (b) and (c) relative to (a). This may 
suggest a preference for either a large interference or a large isovector contribution to the 
subleading terms, or else a mixture of the two. However, in light of the systematic and 
model dependence uncertainties, little can be said with any degree of certainty about the 
nature of the non- diffract ive terms. Leading neutron measurements and data at different 
centre of mass energies would clarify the situation. There are also further uncertainties; 



pion exchange may be relevant in the fixed target data ||2], and measurements at HERA 
with tagged final state protons would be necessary to eliminate systematic effects arising 
from the correction from M^ < 1.6 GeV to the elastic proton vertex. All of the fits give 
a IRIRIR contribution that is consistent with zero, in agreement with previous analyses of 



pp and pp data P2|,|3,H8|,P 



The magnitude of the diffractive component at low W is difficult to constrain in light of the 
uncertainties. In the fits presented here, the diffractive contribution to the integrated cross 
section for 4.0 GeV^ < M^ < 0.05 W^ varies between 1.3 /ib and 2.7 fihatW = 14.3 GeV, 
depending on the assumptions regarding the nature of the subleading exchange. It is 
interesting to note that the value of 3.6 ±0.6 /ib, obtained from a fit with a more empirical 
treatment of non-diffractive contributions performed by [^, is larger than any of the 
values found in this analysis. 

The diffractive contribution at HERA energies and the pomeron intercept are significantly 
better constrained. A single value for the pomeron intercept is obtained by taking an 
average of the three fits. The resuh is a^{0) = 1.068 ± 0.016 ± 0.022 ± 0.041, where the 
errors are statistical, systematic and model related respectively. The spread of values in 
the three different fits is included in the model dependence uncertainty and represents 
the dominant contribution to the quoted error. The value of ^^(O) obtained is similar 
to those extracted from hadronic cross sections that are governed by a soft pomeron 



>(0) ~ 1.081 H). 



In analyses of pp dissociation cross sections extending to centre of mass energies an order 



of magnitude in excess of those investigated here |15|, evidence was found of a need for 



unitarity corrections in the form of a very shallow centre of mass energy dependence. The 
M^ dependence was found to be described by the same pomeron intercept as the total 



cross section. The uncertainties associated with the secondary exchanges in this analysis 
prohibit firm conclusions regarding screening effects for photoproduction. Although fit 
(a), in which neither isovector exchange nor interference is considered, yields a value of 
ajp(O) in the region that might be expected if screening corrections were important ||55|| , 
the value is larger in fits (b) and (c) and in all three fits presented here, the W"^ and M^ 
dependences are well described with a single value for a;jp(0). 
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The diffractive cross section a^ for the process 7p — ^ Xp with 4 GeV^ < M^ < 0.05 W'^ 
is obtained at {W) = 187 GeV by integrating the diffractive terms (PPP and PPIR) 
in equation (|^), with the values for the couphng constants and pomeron intercept taken 
from the fits hsted in table ^. The results of the three fits are averaged to obtain 

a^(4 GeV^ < M^ < 0.05 W^) = 12.4 ± 0.6 (stat.) ± 1.4 (syst.) ± 1.7 (model) /xb . (11) 

Since there is negligible non-diffractive contribution to the cross section in the region 
M^ < 4 GeV^, which is omitted in the fits, the corresponding measured cross section is 
added to that in equation ([TI|) to obtain the integrated cross section for the sum of elastic 
and single photon dissociation processes. The correction from M^ < 1.6 GeV to the fully 
elastic case is again obtained from an average of the two Monte Carlo simulations after 
tuning to match the data. A correction factor of 0.91 ± 0.06 is applied in the second M^ 
interval and 0.96 ± 0.04 in the first. The latter figure is smaller than the correction used 
for the remaining M^ intervals because of the fast fall-off with increasing My in the low 



M^ region, apparent in figure ^ and discussed in section |6.4| . A correction of 1.4 ± 1.4 /ib 



estimated from the di-pion line shape [^ , is added to account for the component of the 
cross section with M^ < 0.4 GeV, which is not measured. The integrated diffractive 
proton-elastic cross section is found to be 

a^{Ml < 0.05 W^) = 34.9 ± 0.9 (stat.) ± 2.3 (syst.) ± 2.7 (model) /ib . (12) 

This result is consistent, within errors, with the previous HI measurement |2^. The diffract 



ive contribution from the process 'jp — >■ Xp to the total photoproduction cross section for 
Ml < 0.05 W^ and (W) = 187 GeV is 22.2 ± 0.6 (stat.) ± 2.6 (syst.) ± 1.7 (model) %. 



6.4 The proton-dissociation cross section 

No previous data or firm theoretical predictions exist that can easily be related to the 
measurement with 1.6 < My < 15.0 GeV. When both M^ and M^ are large a triple-Regge 
analysis is in principle possible @,0, though even with only two trajectories there are a 
total of twelve terms with distinct dependences on M^ , M^ and W"^. The measured double 
dissociation cross section is integrated over a large region in M^ . A full decomposition into 
diffractive and non-diffractive components is not practical without additional assumptions 
as to the ill-constrained behaviour of the cross section in the low M^ region, the relative 
coupling strengths of various reggeons to the photon and the proton and the t dependence 
of triple-Regge couplings. However, as discussed below, certain features of the large M^ 
measurements made here are as might be expected from the conclusions obtained in the 
more detailed analysis of the proton-elastic cross section. 

From the measurements for (W) = 187 GeV with M^ < 1.6 GeV it is found that the 
diffractive cross section in the elastic region is significantly larger than that integrated 
over larger values of M^ (compare equations ( [TI| ) and (|r^) ). A similar pattern is ob- 
served in the dependence on M^ at low M^. The cross section integrated over the range 
0.4 < M^ < 1.26 GeV and M^ < 1.6 GeV is around 3.5 times larger than that for the 
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same M^ region and 1.6 < My < 15.0 GeV. This clearly demonstrates the suppression of 
the cross section for light vector meson production when the proton dissociates compared 
to that when it remains intact. A similar behaviour is observed in the ratios of cross 
sections for the processes pp -^ pp and pp -^ pY |]15|,|4^ . 



Non-diffractive exchanges are likely to become more important with increasing My as 
well as with increasing M^ . This assertion is supported by the observation that the cross 
section for 1.6 < M^ < 15.0 GeV falls less steeply with increasing M^ than that for 
My < 1.6 GeV (see figure ^. In contrast to the vector meson region, the differential 
cross sections in the two ranges of M^ are similar in magnitude for M^ ~ 5 GeV. With 
the present experimental accuracy the large My cross section is well described in the 

measured region by a power law of the form -Trfa- oc ( jj2- ) . In fits to this expression 



dM| - \AP^ 



for 4.0 < Ml < 86.2 GeV^, a value n = 0.84 ± 0.06 (stat.) ± 0.18 (syst.) is obtained at 
(W) = 187 GeV and n = 0.85 ± 0.14 (stat.) ± 0.18 (syst.) at {W) = 231 GeV. 



7 Summary 

A measurement of the cross section da/dM^ for the inclusive process 7p -^ XY has 
been made at a centre of mass energy more than an order of magnitude larger than 
previously. A clear peak in the cross section is observed at the lowest values of M^ and 
My, corresponding to the elastic reaction 7p -^ Vp. The cross section falls rapidly in 
this region as either M^ or M^ is increased. At larger values of M^ an approximate 
dependence dcr/dM^ ~ 1/M^ is observed, both where Y is dominantly a single proton 
and where it has larger masses, though the fall in the cross section is steeper when M^ is 
small. 

For the data with small M^ and moderately large M^, a triple- Regge decomposition of 
the cross section has been presented. Despite the simplifications necessary, it has been 
shown that triple-Regge parameterisations provide a good description of both the W^ and 
the M^ dependence of available photoproduction data, with parameters constrained by 
other processes to which Regge phenomenology is applicable. The data fitted do not lie 
in the asymptotic regime in which only the pomeron need be considered. An important 
feature of the model is that it provides a self-consistent treatment of both diffractive and 
non-diffractive contributions. 

In various model scenarios investigated, the dominant exchange at VT ~ 200 GeV and 
small M^ is found to be diffractive. At the largest dissociation masses accessed, the need 
for a non-diffractive contribution is apparent, which may be described with an effective 
trajectory close to «]j^(t) = 0.55 + 0.90t. More measurements at different centre of mass 
energies or particle identification in the system Y are required if the precise nature of the 
subleading exchanges is to be determined. 

The pomeron intercept is found to be in good agreement with values extracted from total 
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and elastic hadronic cross sections. In global fits to HI and fixed target data, a value 
a^(0) = 1.068 ±0.016 (stat.) ±0.022 (syst.) ±0.041 (model) is obtained. The 'model de- 
pendence' error arises predominantly from the spread of values when different scenarios are 
considered for the subleading terms. The cross section for diffractive processes of the type 
7P -^ Xp for Ml < 0.05 W^ is found to be 34.9 ±0.9 (stat.) ±2.3 (syst.) ±2.7 (model) /ib 
at (W) = 187 GeV, representing 22.2 ± 0.6 (stat.) ± 2.6 (syst.) ± 1.7 (model) % of the 
total photoproduction cross section. 



The double dissociation cross section in the range 1.6 < My < 15.0 GeV and 
4.0 < Ml < 86.2 GeV^ for lU < W < 251 GeV behaves approximately as da/dM^ ~ 

(l/Miy. Yoi{W) = 187 GeV a fit to this form yields n = 0.84±0.06(stat.)±0. 18 (syst.). 
The shallowing of this dependence relative to that for the proton-elastic case may be re- 
lated to large non- diffractive contributions. 
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